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An exact analytical solution of the statistical multifragmentation model is found 
in thermodynamic limit. Excluded volume effects are taken into account in the ther- 
modynamically self-consistent way. The model exhibits a 1-st order phase transition 
of the liquid-gas type. An extension of the model including the Fisher's term is also 
studied. The possibility of the second order phase transition at or above the critical 
point is discussed. The mixed phase region of the phase diagram, where the gas 
of nuclear fragments coexists with the infinite liquid condensate, is unambiguously 
identified. The peculiar thermodynamic properties of the model near the boundary 
between the mixed phase and the pure gaseous phase are studied. The results for 
the caloric curve and specific heat are presented and a physical picture of the nuclear 
liquid-gas phase transition is clarified. 
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Nuclear multifragmentation is one of the most interesting and widely discussed phe- 
nomena in intermediate energy nuclear reactions. The statistical multifragmentation 
model (SMM) (see 0, 0] and references therein) was recently applied to study the liquid- 
gas phase transition in nuclear matter || |], ^ ||. Numerical calculations within the 
canonical ensemble exhibited many intriguing peculiarities of the finite multifragment 
systems. However, the investigation of the system's behavior in the thermodynamic limit 
was still missing. Therefore, there was no rigorous proof of the phase transition existence, 
and the phase diagram structure of the SMM remained unclear. Previous numerical 
studies for the finite nuclear systems (the canonical and microcanonical ensembles) led 
to unjustified (and sometimes wrong) statements concerning the nuclear liquid-gas phase 
transition in the thermodynamic limit. In our recent paper [7f an exact analytical solution 
of the SMM was found within the grand canonical ensemble which naturally allowed to 
study the thermodynamic limit. The self-consistent treatment of the excluded volume 
effects was an important part of this study. In this letter we investigate the peculiar ther- 
modynamic properties near the boundary between the mixed phase and the pure gaseous 
phase. New results for the caloric curve and the specific heat are presented and a phys- 
ical picture of the nuclear liquid-gas phase transition in SMM is clarified. This physical 
picture differs from the one advocated in the previous numerical studies. 

In the SMM the states of the system are specified by the multiplicity sets {n^} (rik = 
0,1,2,...) of fc-nucleon fragments. The partition function of a single fragment with k 
nucleons is [![]: = V (mTk/2ir) 3 ^ 2 Zk , where k = 1,2,..., A (A is the total number 
of nucleons in the system). V and T are, respectively, the volume and the temperature 
of the system, m is the nucleon mass. The first two factors in u>k originate from the 
non-relativistic thermal motion and the last factor, Zk, represents the intrinsic partition 
function of the fc-fragment. For k = 1 (nucleon) we take Z\ — 4 (4 internal spin-isospin 
states) and for fragments with k > 1 we use the expression motivated by the liquid drop 
model (see details in Ref. JTJ): z fc = exp(— fk/T), with the fragment free energy 

fk = ~ [W + T 2 /e Q }k + a{T) k 2/3 + rTlnk. (1) 

Here W Q = 16 MeV is the bulk binding energy per nucleon, T 2 /e Q is the contribution of 
the excited states taken in the Fermi-gas approximation (e Q = 16 MeV) and cr(T) is the 
temperature dependent surface tension which is parameterized in the following form: 

a{T) = a Q [(T 2 - T 2 )/(T 2 + T 2 )f\ (2) 

with (x = 18 MeV and T c = 18 MeV (a = at T > T c ). The last Fisher's term 
in Eq. (0) with dimensionless parameter r is introduced for generality. The canonical 
partition function (CPF) of the ensemble of nuclear fragments has the following form: 

Z^(V,T) = E II -£*"*)■ (3) 

{n k } k=l Uk - k 

The model defined by Eqs.(|T]||) with r = was studied numerically in Refs. || f|, ||, 
[|. This is a simplified version of the SMM since the symmetry-energy and Coulomb 
contributions are neglected. However, its investigation appears to be very important for 
understanding the physics of multifragmentation. 

In Eq. (^) the nuclear fragments are treated as point-like objects. However, these 
fragments have non-zero proper volumes and they should not overlap in the coordinate 



space. In the Van der Waals excluded volume approximation this is achieved by replacing 
the total volume V in Eq. (^) by the free (available) volume Vf = V — bJ2k k n k, where 
b = 1/po (p = 0.16 fm~ 3 is the normal nuclear density). Therefore, the corrected CPF 
becomes: Z A (V, T) = Zf(V - bA, T). 

The calculation of Z A (V,T) is difficult because of the constraint J2k^ n k — A. This 
difficulty can be partly avoided by calculating the grand canonical partition function: 

oo 

Z(V, T,p) = Yl exp (pA/T) Z A (V, T) Q(V - bA) , (4) 

A=0 

where the chemical potential p is introduced. The calculation of Z is still rather difficult. 
The summation over the sets {n^} in Z A cannot be performed analytically because of the 
additional A-dependence in the free volume Vf and the restriction Vf > 0. The problem 
can be solved by introducing the so-called isobaric partition function (IPF) which is 
calculated in a straightforward way (see details in Refs. |7|, |8|, |9|, |i~0fl): 

Z(s,T,p) = / dV exp (sV )Z(V,T,p) = — , (5) 

Jo s — / s.i.u 



where 
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k=2 \ 

with v = p + W Q + T 2 /e . In the thermodynamic limit V — > oo the pressure of the system 
is defined by the farthest-right singularity, s*(T,p), of the IPF Z(s,T,p) 

In Z(V, T, u) ^ . . 

p(T,p) = T im = T s*(T,p) . (7) 

V — ^oo V 

The study of the system's behavior in the thermodynamic limit is therefore reduced to 
the investigation of the singularities of Z. 

The IPF (^]) has two types of singularities: 1) the simple pole singularity defined by 
the following equation s g (T,p) = F(s g ,T,p) ; 2) the singularity of the function T itself 
at the point si(T,p) = vjTb where the coefficient in linear over k terms of the exponent 
in Eq. @ is equal to zero. 

The simple pole singularity corresponds to the gaseous phase where pressure p g = Ts g 
is determined by the following transcendental equation: p g (T,p) = T!F{p g /T : T, p). The 
singularity si(T,p) of the function T defines the liquid pressure: pi(T,p) = Tsi(T,p) = 
vjb. Here the liquid is represented by an infinite fragment (condensate) with k = oo. 

In the region of the (T, /i)-plane where v < bp g (T,p) the gaseous phase is realized 
{.Pg > Pi)i while the liquid phase dominates at v > bp g (T,p). The liquid- gas phase 
transition occurs when the two singularities coincide, i.e. s g (T,p) = si(T,p). As T in 
Eq. (|G|) is a monotonously decreasing function of s the necessary condition for the phase 
transition is that the function JF is finite in its singular point S;. At r = this condition 
requires cr{T) > and, therefore, T < T c . Otherwise, JF(s;,T, p) = oo and the system 
is always in the gaseous phase as s g > sj. As one can see from Eq.(|]) the convergence 



properties of ^(s, T, p) depend significantly on the Fisher's exponent r in the vicinity of 
the critical point where the surface term vanishes. In what follows we mainly concentrate 
on the case r = 0. Other possibilities which appear at r > are discussed shortly. Their 
detail study can be found in Ref. 0. Here we only note that Eqs. ([5], ^ represent an 
exact solution of the Fisher's droplet model fill where additionally the effects of excluded 
volume are incorporated. 

The baryonic density p in the liquid and gaseous phases is given by the following 
formulae, respectively: 

Pi = (9pi/dp) T = 1/6 , p g = (dp g /dp) T = Pid/{1 + bp id ) , (8) 

where the function is the density of point-like nuclear fragments with shifted, p — > 
p — bp g , chemical potential: 



Pu(T,p) = {^)' 



Z\ exp 



' p - bp 
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+ ][> 5 / 2 exp 

k=2 
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(v — bp g )k — ak 2 ^^ 
T 



(9) 



At T < T c the system undergoes a 1-st order phase transition across the line p* = 
p*(T) defined by the condition of coinciding singularities: si = s g , i.e., pi = p g . The phase 
transition line p*{T) in the (T, /x)-plane corresponds to the mixed liquid and gas states. 
This line is transformed into the finite mixed-phase region in the (T, p)-plane shown in 
Fig. 1. The baryonic density in the mixed phase is a superposition of the liquid and gas 
baryonic densities: p = Xpi + (1 — \)p g , where A (0 < A < 1) is the fraction of the system's 
volume occupied by the liquid inside the mixed phase. Similar linear combinations are also 
valid for the entropy density s and the energy density e with (i = l,g) s» = (dpi/dT)^ , 
e % = T (dpi/dT)^ + p (dpi/dp) T - p^ 

Inside the mixed phase at constant density p the parameter A has a specific temper- 
ature dependence shown in Fig. 2: from an approximately constant value p/p at small 
T the function A(T) drops to zero in a narrow vicinity of the boundary separating the 
mixed phase and the pure gaseous phase. This corresponds to a fast change of the config- 
urations from the state which is dominated by one infinite liquid fragment to the gaseous 
multifragment configurations. It happens inside the mixed phase without discontinuities 
in the thermodynamical functions. 

An abrupt decrease of A(T) near this boundary causes a strong increase of the energy 
density as a function of temperature. This is evident from Fig. 3 which shows the caloric 
curves at different baryonic densities. One can clearly see a leveling of temperature at 
energies per nucleon between 10 and 20 MeV. As a consequence this leads to a sharp peak 
in the specific heat per nucleon at constant density, c p (T) = (de/dT) p /p , presented in 
Fig. 4. A finite discontinuity of c p (T) arises at the boundary between the mixed phase 
and the gaseous phase. This finite discontinuity is caused by the fact that A(T) = 0, but 
(d\/dT) p 7^ at this boundary (see Fig. 2). 

The negative values of the specific heat shown in Fig. 4 appear due to the param- 
eterization of the surface tension (FJ): its second derivative with respect to temperature 
generates a negative contribution in the vicinity of T c for all densities. This feature of the 
model is unphysical and it has to be modified. 



It should be emphasized that the energy density is continuous at the boundary of 
the mixed phase and the gaseous phase, hence the sharpness of the peak in c p is entirely 
due to the strong temperature dependence of A(T) near this boundary. Moreover, at any 
p < p the maximum value of c p remains finite and the peak width in c p (T) is nonzero 
in the thermodynamic limit considered in our study. This is in contradiction with the 
expectation of Refs. [|], [5| that an infinite peak of zero width will appear in c p (T) in this 
limit. Also a comment about the so-called "boiling point" is appropriate here. This is a 
discontinuity in the energy density e(T) or, equivalently, a plateau in the temperature as 
a function of the excitation energy. Our analysis shows that this type of behavior indeed 
happens at constant pressure, but not at constant density! This is similar to the usual 
picture of a liquid-gas phase transition. In Refs. [Q, || a rapid increase of the energy 
density as a function of temperature at fixed p near the boundary of the mixed and 
gaseous phases (see Fig. 3) was misinterpreted as a manifestation of the "boiling point" . 

The results presented in Figs. 1-4 are obtained for r = 0. New possibilities appear 
at non-zero values of the parameter r. At0<r<5/2 the qualitative picture remains 
the same as discussed above, although there are some quantitative changes. For r > 5/2 
the condition jF(s,T, p) < oo is also satisfied in the singularity point si(T,fi) for all 
T > T c where <j(T) = 0. Therefore, the liquid-gas phase transition extends now to all 
temperatures. Its properties are, however, different for r > 7/2 and for r < 7/2 (see 
Fig. 5). If t > 7/2 the gas density is always lower than 1/6 as p^ is finite. Therefore, the 
liquid-gas transition at T > T c remains the 1-st order phase transition with discontinuities 
of baryonic density, entropy and energy densities. 

At 5/2 < t < 7/2 the baryonic density of the gas in the mixed phase, p™ %x = 
P?d X ( T )/(. 1 + b p7i x (T)), becomes equal to that of the liquid at T > T c , i.e., p r a ' 
1/6 = p , since 



-mix 
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pT{T) = p id (T,pi*(T)) = (^) 



3/2 



z\ exp — — 



+ £fc!-exp(-^^) ^oo, (10) 

k=2 \ 

if surface tension vanishes a = 0. It is easy to prove that the entropy and energy densities 
for the liquid and gas phases are also equal to each other. There are discontinuities only in 
the derivatives of these densities over T and /z, i.e., p(T, /i) has discontinuities of its second 
derivatives. Therefore, the liquid-gas transition at T > T c for 5/2 < r < 7/2 becomes the 
2-nd order phase transition. According to standard definition, the point T = T c , p = 1/b 
separating the first and second order transitions is the tricritical point. One can see that 
this point is now at a finite pressure. 

It is interesting to note that at r > the mixed phase boundary shown in Fig. 5 is not 
so steep function of T as in the case r = presented in Fig.l. Therefore, the peak in the 
specific heat discussed above becomes less pronounced. 

In conclusion, the simplified version of the SMM is solved analytically in the grand 
canonical ensemble. The progress is achieved by reducing the description of phase tran- 
sitions to the investigation of the isobaric partition function singularities. The model 
clearly demonstrates a 1-st order phase transition of the liquid-gas type. The considered 
system has peculiar properties near the boundary of the mixed and gaseous phases. The 
rapid change of the thermodynamical functions with T at fixed p takes place near this 




boundary due to the disappearance of the infinite liquid fragment. This leads to leveling 
of the caloric curves shown in Fig. 3 at temperatures between 6-10 MeV depending on 
the density. As a consequence a sharp peak and a finite discontinuity are developed in 
the specific heat c p (T) at the boundary of the mixed and gaseous phases. 

The phase diagram appears to be rather sensitive to the value of the parameter r in 
the Fisher's free energy term included in our treatment. New interesting possibilities for 
the phase diagram emerge for r > 5/2 in comparison with the case when r < 5/2. The 
case 5/2 < r < 7/2 is particularly interesting because of the appearance of the tricritical 
point separating the 1-st and 2-nd order phase transitions. 



Acknowledgments. 

The authors are grateful to D. Blaschke, J. Polonyi and P.T. Reuter for stimulating 
discussions. We thank the Alexander von Humboldt Foundation and DAAD (Germany) 
for the financial support. The research described in this publication was made possible in 
part by Award No. UP1-2119 of the U.S. Civilian Research & Development Foundation for 
the Independent States of the Former Soviet Union (CRDF). K.A.B. gratefully acknowl- 
edges the warm hospitality of the Particle- and Astrophysics Group of the University of 
Rostock, where this talk was given. 

References 

[1] J. P. Bondorf, A.S. Botvina, A.S. Iljinov, I.N. Mishustin, K.S. Sneppen, Phys. Rep. 
257 (1995) 131. 



[2 
[3 

[5 
[6 

[7; 



[9 
[10; 

[11 



D.H.E. Gross, Phys. Rep. 279 (1997) 119. 

K.C. Chase and A.Z. Mekjian, Phys. Rev. C 52 (1995) R2339. 

S. Das Gupta and A.Z. Mekjian, Phys. Rev. C 57 (1998) 1361. 

S. Das Gupta, A. Majumder, S. Pratt and A. Mekjian, |nucl-th/9903007| (1999). 



S. Das Gupta, A.Z. Mekjian and M.B. Tsang, |nucl-th/ 0009033] (2000). 

K.A. Bugaev, M.I. Gorenstein, I.N. Mishustin and W. Greiner, Phys. Rev. C62 
(2000) 044320. 

M.I. Gorenstein, V.K. Petrov and G.M. Zinovjev, Phys. Lett. B 106 (1981) 327; 
M.I. Gorenstein, W. Greiner and S.N. Yang, J. Phys. G 24 (1998) 725. 



K.A. Bugaev, M.I. Gorenstein, I.N. Mishustin and W. Greiner, [nucl-th/0007062 

(2000) 3 p. 



K.A. Bugaev, M.I. Gorenstein, I.N. Mishustin and W. Greiner, |nucl-th/00 11055 
(2000) 7 p. 

M.E. Fisher, Physics 3 (1967) 255. 




5 10 15 T c 20 

T (MeV) 



Figure 1: Phase diagram in the (T, p)-plane. The mixed phase and pure gaseous phase 
boundary is shown by the solid line. The pure liquid phase (shown by crosses) corresponds 
to the fixed density p = p Q . Point C is the critical point, at T > T c only the pure gaseous 
phase exists. 
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Figure 2: Volume fraction A(T) of the liquid inside the mixed phase is shown as a function 
of temperature for fixed nucleon densities p/p Q = 1/6, 1/3, 1/2,2/3,5/6 (from bottom to 
top). 




Figure 3: Temperature as a function of energy density per nucleon (caloric curve) is 
shown for fixed nucleon densities p/ p Q = 1/6, 1/3, 1/2,2/3. 




Figure 4: Specific heat per nucleon as a function of temperature at fixed nucleon density 
p/p = 1/3. The dashed line shows the finite discontinuity of c p {T) at the boundary of the 
mixed and gaseous phases. The negative specific heat appears in the vicinity of T = 18 
MeV. 
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Figure 5: Phase diagrams in T — p plane for r = 3.6 (upper panel) and r = 2.6 (lower 
panel). Point C in the lower panel is the tricritical point. Crosses correspond to the liquid 
phase of the first order phase transition and dots correspond to the states of the second 
order one. 



